Anharmonicity induced resonances for ultracold atoms and their detection 
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When two atoms interact in the presence of an anharmonic potential, such as an optical lattice, 
the center of mass motion cannot be separated from the relative motion. In addition to generating 
a confinement-induced resonance (or shifting the position of an existing Feshbach resonance), the 
external potential changes the resonance picture qualitatively by introducing new resonances where 
molecular excited center of mass states cross the scattering threshold. We demonstrate the existence 
of these resonances, give their quantitative characterization in an optical superlattice, and propose 
an experimental scheme to detect them through controlled sweeping of the magnetic field. 
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I. INTRODUCTION 

In recent years, there has been much progress in the 
study of ultracold atoms in optical lattices, which can 
cleanly emulate important models in condensed matter, 
hold promise for quantum computing schemes, and offer 
the prospect to observe many interesting new phenom- 
ena The versatility of this line of research is due in 
no small part to the control of the atomic interactions 
afforded by tuningan external magnetic field near a Fes- 
hbach resonance [2|. In addition to a magnetic field, a 
confining potential can also be used to tune the scat- 
tering length via a Feshbach-type mechanism, typically 
referred to as a confinement-induced resonance [3[ or a 
trap-induced shape resonance [|| depending on the trap 
configuration. The trap-induced resonance is basically 
caused by a shift of the free-space Feshbach resonance 
point by the confining potential In an optical lattice, 
the possibility of decay of atomic pairs due to anharmonic 
coupling to molecules in an excited center-of-mass (cm.) 
state has previously been discussed in Ref. [f|. More re- 
cently, the possibility of a controlled transfer of an atomic 
pair to a molecule in an excited cm. state was discussed 
in Ref. Q ■ The anharmonicity of the optical lattice po- 
tential has also been recognized as important in obtain- 
ing quantitatively accurate predictions for the shift of the 
free space Feshbach resonance position, binding energy, 
etc. [|. 

In this paper we point out a new effect whereby an- 
harmonic confinement, e.g., from an optical lattice, not 
only shifts the free-space resonance point, but also in- 
duces a series of additional scattering resonances. (A 
similar effect occurs in mixed dimensions in the absence 
of anharmonicity 0.) Thus, anharmonicity may give 
rise not only to population transfer between atom pairs 
and molecules different cm. states, as mentioned to var- 
ious extents in previous works [6|— Isj , but also to a strong 
modification of the effective atom-atom interaction in the 
vicinity of the induced scattering resonances. In order to 
use the optical lattice system as a quantum emulator, it is 
important to have a full understanding of the dependence 
of the interaction on the experimental parameters. Even 
such a basic item as the form of the effective many-body 



(a) 




molecular 
continuum 



atomic 
continuum 



free space 
Feshbach resonance 



(b) 4 




induced 
resonance 



f~ molecular 
vv^J))) excited band 



molecular 
ground band 

atomic 
ground band 



shifted free 
space resonance 



(c) 



,((#)) 




induced 
avoided crossing 



molecular 
excited state 
molecular 
ground state 

atomic 
ground state 



shifted 
avoided crossing 



magnetic field 



FIG. 1: Sketches of the Feshbach type of resonances (a) in free 
space; (b) in an optical lattice with additional anharmonicity 
induced resonances; (c) in a confining potential where the 
resonances are signaled by the avoided level crossings. 



lattice Hamiltonian [Tl| will be affected in the vicinity 
of an induced resonance (for a detailed treatment, see 
Ref. [I3|). The presence of the additional resonances is 
then an important consideration for experiment as well 
as a novel tool for tuning the interaction utilizing a res- 
onance between atoms and excited cm. molecules. Mea- 
suring population transfer between cm. states can pro- 
vide a handy means to look for the resonances. Below, 
we characterize these anharmonicity induced resonances 
in an optical superlattice and propose an experimental 
scheme to detect their consequences. 

To understand the basic mechanism of the anhar- 
monicity induced resonances, let us first compare it with 
the free-space Feshbach resonance. The free space Fesh- 
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bach resonance is caused by coupling between the scatter- 
ing state of the atomic pair and a highly excited molecu- 
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lar level (the Feshbach molecule), as depicted in Fig. 1(a) 
When the energy of the Feshbach molecule, tuned by 
the external magnetic field, crosses the lowest scattering 
state, a resonance in the scattering length is signaled ■ 
In free space, cm. and relative motions are decoupled 
during the atomic scattering, and the cm. momentum 
forms a continuum which is not altered by the scattering 
process. 

In the presence of an optical lattice, the continuum 
spectrum for the atomic and the molecular cm. motion 
both split in a series of energy bands. We consider scat- 
tering of the atoms in the lowest bands, and only this 
lowest atomic band is shown in Fig. 1(b) However, even 
for this lowest-band atomic scattering, the excited bands 
for the cm. motion of the Feshbach molecule still play a 
significant role due to the anharmonicity of the optical 
lattice potential. In a harmonic potential, the cm. mo- 
tion of two colliding atoms is separated from their relative 
motion, and thus remains in the lowest band during the 
collision and does not couple to the Feshbach molecule 
in the excited bands. However, the anharmonicity of the 
potential mixes the cm. and relative motions, and the 
lowest band scattering state of the atoms is coupled to the 



Feshbach molecule in each band, as depicted in Fig. 1(b) 
As one can see from this figure, all the bands for the Fes 
hbach molecule, no matter how excited, eventually cross 
the atomic scattering threshold as the magnetic field is 
lowered. This will lead to many resonances for the atomic 
scattering. In practice, the anharmonic coupling between 
a Feshbach molecule in the excited band and the atomic 
pair state in the lowest band will decrease as the band 
becomes more excited, and the resonances become pro- 
gressively narrower as the magnetic field is lowered, so 
only the first few of these resonances are broad enough 
to be experimentally observable. 

In order to quantitatively characterize the anharmonic- 
ity induced resonances, we consider the atoms in an op- 
tical supcrlatticc potential. In an optical lattice, direct 
calculation of the scattering length between two atoms 
is challenging as one can not separate the cm. and the 
relative motion and solution of an equation with all six 
degrees of freedom is numerically demanding. Instead, 
here we consider the atoms in a deep superlattice poten- 
tial [lH , which separates the periodic optical lattice into 
a series of double-well potentials. This has several moti- 
vations: First, by adding a confining trap, as illustrated 
in Fig 
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the resonance in the continuum scattering 
spectrum caused by the emergence of a new Feshbach 
molecular level becomes an avoided level crossing in the 
discrete spectrum of the trapped atoms. By calculating 
the width and the position of the avoided level crossing, 
we can approximately characterize the resonance prop- 
erties for the atomic scattering. Numerically, it is more 
convenient to deal with the discrete spectrum in a trap 
which allows application of specific calculation techniques 
presented below. Second, the optical superlattice potcn- 
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FIG. 2: (Color online.) (a) Double-well potential (solid) for 
a single atom along z and an example of a corresponding 
superlattice potential (dashed). The horizontal lines are the 
lowest four single-atom energy levels of the double-well, (b) 
Contour plot of the locally isotropic 3D double- well potential. 



tial has been realized in experiments [l2j , which allows 
direct detection of consequences of the anharmonicity in- 
duced resonances in this kind of trap. We will propose an 
experimental scheme to test the quantitative predictions 
from the anharmonicity induced resonances in a superlat- 
tice. Third, the anharmonicity induced resonances also 
affect the effective many-body Hamiltonian for strongly 
interacting atoms in an optical lattice [l(J [H| • A natural 
step to derive such a Hamiltonian is to first consider the 
effective interaction for atoms in double-well potentials 
realized with a deep optical superlattice. 



II. METHODS 

We assume that the superlattice potential is along the 
axial direction z which separates the system into a se- 
ries of double- wells [HI]. We consider two atoms of mass 
to in each double- well potential V(z), approximated by 
Taylor expanding Vq cos 2 {Ulz) to 12 th order in z. Here 
Vq sets the barrier depth and k^ sets the distance be- 
tween wells, and is related to the laser wavevector. Al- 
though it is not important for our purposes to exactly fit 
a particular form of superlattice potential, if one takes 
a superlattice of the form cos 2 (kz) 
one should choose kh - 



csin (kz/2), then 
as shown in 
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Fig. [2j (We express energy in units of a "recoil energy," 
Er = /i 2 fc 2 /2m, and plot the case Vb = &Er.) In any 
case, this potential should be quite sufficient to capture 
the essential physics in the limit of independent double- 
wells. For ease of calculation, the lattice wells in the 
transverse directions are approximated by harmonic po- 
tentials, with the frequency, oj, chosen such that the po- 
tential is locally isotropic at the bottom of each well. 

Due to the harmonicity of the transverse trap, the cm. 
motion in the transverse direction separates out and is 
thus neglected in the rest of the discussion. Also, due to 
the axial symmetry of the trap, the azimuthal angular 
dependence of the relative motion separates out. How- 
ever, along the axis of the double-well, the cm. motion 
is not separable from the relative motion. The two atom 
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system then has three relevant coordinates: the relative 
coordinates z and p, along the axial and transverse di- 
rections, respectively, and the axial cm. coordinate. 
In terms of these coordinates, the external potential with 
barrier depth V is V(p, z, Z) = V(p) + V(z, Z), where 

V(p) = V k 2 L p 2 /2 

v(z z)-vY y (~4)"r(i-2n) 2w / z yn 

V(z, Z) - V ^ r (1 _ 4n) r (1 + in) k L \Z±-J 

with r (x) the Eulcr gamma function, and our summation 
over signs denotes that for each value of n one must also 
add the two terms corresponding to the upper and lower 
signs. 

The atoms are interacting via a short-range poten- 
tial U(r) characterized by its s-wave scattering length 
a s . The exact form of the interaction is irrelevant as 
long as its effective range is much smaller than the aver- 
age interatomic distance and the trap length scale. For 
a broad s-wave Feshbach resonance, the use of a zero- 
range pseudopotential is typically justified |13T ]. as has 
also been confirmed experimentally [Tij . Numerically, it 
is easier to use a finite-range attractive Gaussian interac- 
tion U(r) = — Uq exp (— r 2 /Vg) , where we typically take 

ro = Q.0'dy , h/muj. Finite-range effects should be neg- 
ligible for such small values of tq, and we have verified 
this by repeating our calculations with tq = O.l^h/muj. 
The free space scattering length is varied by adjusting the 
strength of the interaction, Uq. We have used values of Uq 
such that the potential supports either zero bound states 
(for negative scattering length) or one bound state (for 
positive scattering length). The scattering length goes 
through resonance when the lowest eigenstate of the in- 
teraction potential passes from being unbound to bound. 

Adopting units such that fc^ = 1 and Er = 
K 2 k\/2m = 1, the Hamiltonian may be written as 

2 d d d 2 Id 
p ^ 

p dp dp dz 2 2 dZ 2 

+ 2 1 4 + V(p,z,Z)-U Q e-(^y r o (1) 
P 

where rri£ is the relative angular momentum, which is a 
good quantum number due to axial symmetry. In the 
following we will only consider mt = 0, since, in the limit 
as ro goes to zero, the interaction does not affect states 
with rrif ^ 0. 

We find the low-lying states of the system using a 
stochastic variational method [l5[ . In this approach, the 
variational wavefunction takes the form 

JV 

*(p, z,Z) = Y / on exp(-p 2 /a 2 - z 2 /b 2 - Z 2 /c 2 ) , (2) 

i 

where a is a linear variational parameter, {a, b, c} are 
nonlinear variational parameters which define the basis 



elements, and N is the size of the basis set. The non- 
linear parameters are selected from stochastically gener- 
ated pools of candidates to minimize the variational en- 
ergy | i/ 1 The basic algorithm is as follows: 
starting with a set of N — 1 basis states, 

1) a pool of (in our calculations) 25 new basis states is 

randomly generated, each defined by a given value 
of {a,i, bi,Ci}; 

2) for each of the 25 possible ./V-dimensional basis sets 

formed by adding one basis from the candidate 
pool, the energy is minimized with respect to a; 

3) the new basis set that yields the lowest energy is kept 

and the previous steps are repeated until the basis 
size, N, increases to the desired number. 

Once every few iterations, the existing basis set is op- 
timized by the following refining process: starting with a 
set of N basis states and n = 1, 

A) a pool of 25 replacement basis states is randomly gen- 

erated, each defined by a given value of {a„, b n ,c n }; 

B) for each of the 25 possible A^-dimensional basis sets 

formed by replacing the n th old basis state with 
a new one from the candidate pool, the energy is 
minimized with respect to a; 

C) if the lowest of these 25 energies is lower than the 

current variational energy, the n th old basis state 
is replaced by the new optimal one and the previous 
steps are repeated for n = 1...N. 

We typically achieve fairly good convergence for N ~ 
300. Although in principle the nonlinear basis optimiza- 
tion must be performed for each value of a s , actually 
the basis set does not change too much as one sweeps 
across resonance except to include narrower and narrower 
Gaussians for positive a s where deeply bound molecules 
form. Apart from deeply bound states, the change in 
the wavefunction is mainly due to changing the expan- 
sion coefficients, a. To save computational time then, 
we performed the nonlinear basis optimization for four 
different values of a s ranging from positive to negative, 
joined the four optimized basis sets, and simply mini- 
mized the energy with respect to a using the resultant 
basis set of about 1200 elements for all values of a s . 

III. RESULTS 

A. Spectrum 

In Fig. [31 we show the energy spectrum of two par- 
ticles interacting near a free space Feshbach resonance 
(l/kLCL s — 0) in the double- well potential. For clarity, 
we have omitted the levels corresponding to wavefunc- 
tions of odd parity in z or Z (which have no contribution 
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(a) 




(b) 



FIG. 3: (Color online.) (a) Spectrum of two strongly in- 
teracting atoms in a three-dimensional double-well potential 
with Vo = 6Er. Only the first few plunging levels are shown, 
(b) Close-up of the strongly interacting region. The first few 
states are explicitly labeled for reference. 



to the anharmonicity induced resonances) and plunging 
levels for —l/ki,a s < —10. We have explicitly labeled 
the lowest few states for later reference. To understand 
Fig. [3j it is useful to use the language of the two-channel 
picture of atom pairs coupled to molecules, as in Fig. [T] 
Without coupling, there are plunging molecular levels 
crossing flat (i.e., noninteracting) atom pair levels, as 
depicted by the dashed lines in Fig. 



1(c) When one 



turns on atom-molecule coupling only between molecules 
and atoms with the same cm. motion, the crossings be- 
tween the lowest plunging molecular level and the flat 
atomic levels become avoided crossings and the spectrum 
is similar to the well-known results for a harmonic trap 
[lfjj . The atoms and the lowest cm. molecules hybridize, 
such that as the inverse scattering length is adiabatically 
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swept from negative values to positive values, the lowest 
atomic level evolves into the lowest molecular level and a 
given excited atomic level will evolve into the next low- 
est atomic level, sweeping out a sigmoidal path. (For a 
double- well, the lowest two atomic levels form a closely 
spaced doublet, so in the presence of coupling the lowest 
sigmoidal level is essentially flat, as in Fig. 
lying doublets behave likewise.) 

If we take anharmonicity into account by allowing cou- 
pling also between the atoms and excited cm. molecules, 
the crossings between the higher plunging molecular lev- 
els and the flat atomic levels also become narrow avoided 
crossings. These signal the presence of a rich set of 
induced resonances. The resonances are weak relative 
to the free space resonance, and become progressively 
weaker away from l/fc^a s = 0, so that only the first few 
are observable. Diabatically, then, Fig. [3] displays three 
kinds of curves: plunging levels corresponding to tightly 
bound molecules, flat levels corresponding to atoms in 
separate wells, and sigmoidal levels corresponding to 
atoms with overlapping wavefunctions such that they in- 
teract while maintaining a nonvanishing pair size unlike 
the tightly bound molecules. Note that the many plung- 
ing molecular levels, of which we have shown only the first 
few, are associated with the various states of the trap, as 
sketched in Fig. 



1(c) 



They are motionally excited cm. 
states, not internally excited states of the interaction po- 
tential. Also note that for an optical lattice, one will 
obtain a similar spectrum except that the discrete levels 
of the double- well shown in Fig. [3] (analogous to Fig. 1(c) ) 



will broaden into bands (analogous to Fig. 1 (b) ) 



B. Avoided crossing data 

To characterize the anharmonicity induced resonance, 
we estimate the time required to adiabatically sweep 
across the avoided level crossing, transferring population 
between atomic and molecular states. In the Landau- 
Zcner approximation [17| . the probability of an adi- 
abatic transfer at sweep rate dB/dt = v is P a d = 
1 — cxp (— vlz/v), where the Landau-Zcncr parameter 
vlz = 7rA 2 /2fi,|9A/9-B|, A is the minimum energy gap 
between the two levels in question, and <9A jdB is the rate 
at which the energy gap changes with the magnetic field 
away from the avoided crossing. The energy splitting A 
for the avoided level crossing should be proportional to 
the width of the corresponding anharmonicity induced 
resonance in a periodic optical lattice. This parameter 
A is listed in Table U for the various avoided crossings 
between excited molecular states and the lowest atomic 
level (near E = 9.3Er) for 6 Li or 40 K atoms. The num- 
bers quoted are of course only a rough guide to what may 
be expected in experiment, and are not intended to be 
quantitatively precise - recall that the double-well po- 
tential we have taken is only an approximation to what- 
ever form the actual potential may take. To connect our 
results to experiment, we assume the scattering length 
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— l/khda 


A/ft (kHz) 


tmin (ms) 


VLZ (G/s) 




0.2 


8 


40 


2 x 10 6 




-0.9 


8 


70 


3 x 10 s 


ft T . 

b Li 


-1.4 


5 


100 


6 x 10 4 




-1.9 


1 


200 


2 x 10 3 




-2.5 


0.3 


1 x 10 3 


40 




0.2 


1 


600 


800 




-0.9 


1 


600 


400 


40 K 


-1.4 


0.8 


1 x 10 3 


100 




-1.9 


0.2 


2 x 10 3 


6 




-2.5 


0.05 


1 x 10 4 


0.1 



TABLE I: Anharmonicity induced avoided level crossing data 
for 6 Li ( 40 K) atoms at Vo = 6Er. These are the avoided 
crossings near E = 9.3E R shown in Fig. U 




FIG. 4: (a) Energy gap and (b) Landau-Zener parameter 
for the first four avoided crossings vs. well depth for 40 K 
atoms. From the top to the bottom, the cu rves correspond to 
the avoided crossings shown in Fig. 3(b) at E = 9.3-Er and 
— 1/kLds = 0.2, —0.9, —1.4, —1.9, respectively. 



is related to the magnetic field via the usual relation 
«s = o-bg [1 — W/ (B — Bq)], where db g is the background 
scattering length, W is the resonance width, and Bq is 
the resonance point. We take ~ 2tt /l\xm and consider 
6 Li[^°K) near the free space Feshbach resonance at 834 
G [ill (202 G [lj|). In TableU we have also listed an esti- 
mate of the minimum time, t m i n , required to ramp across 
the avoided crossing at the critical rate, vlz- If the time 
available in the experiment to perform the ramp is on the 
order of a few milliseconds, appreciable adiabatic trans- 
fer is feasible across the first five (four) avoided crossings 
for 6 Li ( 40 K) atoms. 

We have performed the same kind of calculations for 
several lattice depths. In Fig. |||we show how the energy 
splitting A and the Landau-Zener parameter vlz for the 
first few resonances listed in TableU change as Vo is varied 
for 40 K. Generally, the energy splitting for the avoided 
level crossing decreases for deeper wells, as one would ex- 
pect due to suppression of the anharmonicity in a deep 
lattice (the harmonic approximation becomes better for 
a deep lattice well). For very shallow wells, though, the 
potential apparently cannot couple the higher cm. states 
of Feshbach molecules to the lowest atomic state as ef- 



ficiently, and the energy splitting A actually increases 
with the lattice depth at first for small Vo- However, 
for very excited cm. states (corresponding, e.g. to the 
bottom curve in Fig. 2]), the weak coupling to the lowest 
atomic state evidently does not depend as strongly on 
the lattice depth. As the potential wells are deepened, 
the resonance positions shift slightly to lower magnetic 
fields. 



C. Detection 

To experimentally detect the avoided level crossings 
associated with the anharmonicity induced resonances, 
one can take the following steps: First, one loads the op- 
tical supcrlattice in the weakly interacting region with 
two atoms in each double well fl2l . [2?jj . The inter- well 
barrier is kept high so that one has a Mott state with 
one atom per well. Second, one ramps the system to the 
strongly interacting region with — l/fcj,a s = ±2, and then 
quickly lowers the inter-well barrier to the desired value 
(with Vq — QEji in our example), leaving the atoms still 
in the Mott state (at energy E ~ 9.3-Er in Fig. [3]) at 
this moment. The magnetic field is then adiabatically 
ramped across the anharmonicity induced resonances, 
and one detects the resulting population distribution af- 
ter the ramp. To do the detection, the inter-well barrier 
is quickly turned back up with a time scale fast com- 
pared with the inter-well dynamics, but still slow com- 
pared with the lattice band gap (or the single-well en- 
ergy gap). This freezes the system evolution again be- 
fore the magnetic field is ramped to the deep BEC side 
(—1/kLds <C —1), separating the molecular levels from 
the atomic levels. One can then selectively take absorp- 
tion images of either the atoms or the molecules [2l[ , and 
measure their distribution over different bands through 
a band-mapping procedure [l2| . The presence of the an- 
harmonicity induced level crossings can then be inferred 
from the final population distribution. 



As an example, in Fig. 5(a) we show the Landau-Zener 
calculation results for 40 K atoms at Vo — 6Er swept from 
— 1/kLds = —2 to —1/kLds = 2 with the atoms starting 
in the state 6 labeled in Fig. [3(b)] For fast sweeps, the 
atoms change states diabatically to remain at about the 
same energy, as would be expected in the absence of an- 
harmonicity. In the adiabatic limit, all the atoms remain 
in state 6, which corresponds at —l/kLd s = 2 to atoms in 
an excited state. At intermediate speeds, several atomic 
states become populated. A sweep in the opposite direc- 
tion, from —1/kLds = 2 to — l/kLd s = —2, starting with 



atoms in state 2, is shown in Fig. 5(b) When sweeping 
in this direction, population can be transferred to tightly 
bound molecules in several excited cm. states (states 2-5 
at —1/kLds = —2) as well as diabatically to atoms near 
the initial energy (state 6). 
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FIG. 5: (Color online.) Final population distribution 
vs. ramp speed of the magnetic field (a) from state 6 at 
— l/fcz,a s = —2 to states 2-7 at —l/kLa 3 = 2 or (b) from 
state 2 at — 1/feLffls = 2 to states 2-6 at — l/kha s — —2. Both 
plots are for 40 K atoms with Vb = 6-Er. In (b), the proba- 
bility of sweeping into state 5 at — l/kLa 3 = —2 is essentially 
zero at any ramp speed. 



IV. SUMMARY 



We predict the existence of several Feshbach-type res- 
onances induced by the anharmonicity of the optical lat- 
tice, which couples the Fcshbach molecules in the excited 
bands and the atomic states in the lowest band. We have 
characterized the corresponding set of avoided level cross- 
ings in the calculated spectrum of two atoms interacting 
in a superlattice potential, and proposed an experimental 
scheme to observe these avoided crossings through slow 
sweeps of the magnetic field. The anharmonicity induced 
resonances may prove to be a useful tool for manipulation 
of interaction between ultracold atoms in optical lattice 
potentials. 

This work was supported by the AFOSR through 
MURI, the DARPA, and the IARPA. 
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